Spin orbit effects in a GaAs quantum dot in a parallel magnetic field 
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We analyze the effects of spin-orbit coupling on fluctuations of the conductance of a quantum dot 
fabricated in a GaAs heterostructure. We argue that spin-orbit effects may become important in 
the presence of a large parallel magnetic field Bu , even if they are negligble for Bu = 0. This should 
be manifest in the level repulsion of a closed dot, and in reduced conductance fluctuations in dots 
with a small number of open channels in each lead, for large Bu . Our picture is consistent with the 
experimental observations of Folk et al. 

PACS numbers 



Recent experiments by Folk et al. studied statistics 
of fluctuations of the conductance g through a quantum 
dot in a GaAs heterostructure with an applied magnetic 
field B\\ in the plane of the sample. In the largest dots 
studied, the application of B\\ was observed to reduce the 
variance of the fluctuations, var(g), by a factor of roughly 
four, in contrast to a reduction factor of two, which was 
originaly expected. As noted by Folk et al., the extra 
reduction might be understood if, for some reason, spin- 
orbit coupling increased with the application of Bu . 

In this note we argue that there should indeed be field- 
dependent spin-orbit effects which are unique to a quan- 
tum dot, arising from the Rashba and/or Dresselhaus 
terms in the effective hamiltonian for electrons in a GaAs 
heterostructure. These effects could well explain the ob- 
servations of Rcf . 0] . Similar effects should appear in the 
repulsion between energy levels in a closed dot. 

We first explain how spin-orbit coupling and parallel 
magnetic fields affect conductance fluctuations. We then 
show that in quantum dots, unlike open systems, the ef- 
fectiveness of spin-orbit interaction grows with increasing 
_B| | . Finally, we examine the consequences of this effect 
on conductance fluctuations in quantum dots, and its re- 
lation to the measurements in Ref. (!]]. 

In a single particle picture, conductance fluctuations 
through a chaotic or disordered quantum dot may be 
crudely understood as arising from fluctuations in the 
number of electronic levels in an energy window of size 
2NA, and in the matrix elements coupling these levels 
to the leads. Here A is the mean level spacing in the dot, 
for each spin state, and N is the number of channels in 
each lead, (i.e., each lead has conductance 2Ne 2 /h). We 
assume the leads to be perfectly coupled to the dot, such 
that Coulomb blockade effects are insiginificant. The 
mean conductance in this geometry, including both spin 
states, is (g) — Ne 2 /h. In the experiments of Ref. Q|, N 
was in the range 1 to 3. 

In the experiments of Ref. a weak perpendicular 
magnetic field B± was applied. This field was strong 
enough to break time-reversal symmetry for the orbital 
motion, but not strong enough to produce a significant 
Zeeman splitting. Then, if spin-orbit coupling is absent 



and Bu = 0, conductance fluctuations should satisfy 

var(g) = 4Cjv, (no spin-orbit, Bu =0) (1) 

where the constant C/v is var(g) for spinless electrons in 
a dot with N open channels per lead, and the factor 4 
results from the degeneracy of the two spin states. (From 
here on we measure all conductances in units of e /h, 
so that Cjv is dimensionless.) The factor Cjv depends 
on the temperature T through the ratio T/Th (where 
r = NA/nh is the escape rate from the dot) and on 
the phase breaking rate r^" 1 through the parameter Tt^. 
The value of Cjv can be calculated from random matrix 
theory, using the Gaussian Unitary Ensemble (GUE). 

For B\\ 7^0, still in the absence of spin-orbit coupling, 
the Fermi levels for spin-up and spin-down electrons are 
split by the Zeeman energy Ez = g* ij,bB\\. When Ez 
is larger than both T and V, the contributions from the 
two spin states become statistically independent, giving 

var (g) = 2C N (no spin-orbit, E z > I\ T). (2) 

However, in the presence of a strong spin-orbit cou- 
pling, the two spin levels will be mixed, and will be de- 
scribed by a single GUE, with mean level-spacing A/2, 
and 2N open channels in each lead. (Recall that Kramers 
degeneracy is already broken by B±.) Thus in that case, 



var (g) = Q 



2JY 



(strong spin-orbit ). (3) 



The cross-over to strong spin-orbit coupling should be 
controlled by the dimensionless parameter A = e so /A 
where e so is the root-mean-square (RMS) value of the 
matrix clement (i \ H so \j) . Here the states i,j have op- 
posite spin directions and orbital energies that differ by 
Ez- (The matrix element is to be calculated with the 
dot isolated from the leads. The same parameter A con- 
trols the repulsion between levels of opposite spins in the 
closed dot.) Then in the presence of Bu we can write 



var {g) Bu = F N (X,T '/hi \Tt^ 



(4) 



where Fjv — > 2Cjv for A — > 0, and Fn — > C-jjv for A suf- 
ficiently large. Note that r is unchanged if N is doubled 



1 



and A is halved, so T remains constant as one varies A. 
We shall also see that at least approximately, Cm ~ C2N, 
so that Fjv decreases by a factor of two as A varies from 
to oo. Then, if the system parameters are such that 
A grows from zero to a large value as a parallel field Bu 
is turned on, the factor-of-two reduction in F, combined 
with the the factor-of-two reduction on breaking the spin 
degeneracy, should lead to overall reduction of a factor 
of m 4 in var(<?), relative to the B\\ = value, in Eq. (Q). 
This is in accord with the observation of Ref . [EJ . 

A variety of evidence, based on Random Matrix The- 
ory and other approaches, suggests that for a large N, 
Cjv is independent of N, for any fixed value of T, T, and 
T0 The biggest deviation from this is presumably 

for N = 1 and no dephasing. At T = 0, with no dephas- 
ing, the value of Cn is known, within Random Matrix 
Theory, to be (16 - AN- 2 )- 1 §. Thus for N = 1, the 
reduction factor Fjv(0, 0, oo)/-Fjv(oo, 0, oo) is 5/2 rather 
than 2. 

We define a crossover value A c where var(g) is half- 
way between the values for A = and A large. We may 
estimate A c as the value of A such that r^ 1 = t^ 1 + T, 
where r^ 1 = 2ttX 2 A/H is the rate for spin flip scattering 
due to the spin-orbit coupling, given by Fermi's Golden 
Rule. Writing N eS = N + tt^At^)' 1 , this gives 



A c « 0.237V 



1/2 



elf 



(5) 



Numerical calculations, discussed further below and il- 
lustrated in Fig. [1], are at least qualitatively consistent 
with this estimate, but suggest that the factor 0.23 should 
be replaced by 0.1. 

Spin-orbit coupling in GaAs heterostructures origi- 
nates from the asymmetry of the potential creating the 
2DEG (Rashba term) and from the lack of inversion sym- 
metry in the GaAs lattice structure (Dresselhaus term). 
The operator describing the spin-orbit coupling is com- 
posed of both terms: 



: 7(« x a) ■ z + T](v x a x - v y a y ), 



(6) 



where v is the velocity operator, a are the Pauli spin 
matrices, 7 and r] are coupling constants, and we ignore 
terms oc v 3 . We assume that the 2DEG is grown on a 
[001] GaAs plane and x, y denote the cubic axes in the 
plane. In an open clean system, Eq. ([|) leads to a small 
spin- splitting of the conduction bands. Note that this 
spin-orbit coupling is different from the one encountered 
in measurements of conductance fluctuations in metals. 
The latter is induced by impurities, and is characterized 
by coupling constants that strongly vary with position. 

In the absence of H so , the eigenstates of the electronic 
Hamiltonian are products of a spatial part \a) and a spin 
part |cr), where a =f (|) denotes spin parallel (antiparal- 
lel) to the Zeeman field B\\. With H so , the mean-square 
value of the dimensionless spin-orbit coupling A for states 
with opposite spins at the Fermi energy is 



\ 2 (E z )^J2\( H ^)atbl\ S(e a -6 h -E z )5(e a ~e F ) (7) 

ab 

where e a , f-b are the orbital energies of the states \a) and 
\b) (i.e., the energies at Bu =0), the over-bar denotes 
averaging over disorder, and e F is the Fermi energy. Here 
and henceforth O s i jS 2 = (si \0\ s2). As we now show, 
in a quantum dot the typical matrix element (H so ) a -\^i 
depends on the energy difference e a — £b, so that A 2 does 
indeed depend on Ez- 

For simplicity, we first discuss the case where r\ = in 
H SOl and we choose S11 || x. For a macroscopic system 
in the diffusive regime, comparing the Kubo-Greenwood 
formula with the Drude formula, one finds 



\{v x )ab\' 



2DA 



1 



nh 1 + (lot) 2 



(8) 



where ftw = e a — tb, t is the transport lifetime, and D = 
v f t/2 is the diffusion constant. Thus, for a large diffusive 

1 /2 

system and for lot <C 1, we find A w {~k) 1 which 
does not depend on to. A parallel field does not affect the 
strength of spin-orbit coupling as long as Ez <C H/t. 

In contrast, the confinement of the electron to a quan- 
tum dot suppresses the velocity matrix elements when 
lotji <C 1, where tr is a Thouless time, which we define 
as the time for an electron to cross from the center to 
the edge of the dot. This is most easily seen if we use 
the relation |(ur)a&| = w|x Q t,|, and note that the matrix 
element of x is bounded by the maximum radius R of the 
dot. More precisely, we may use the relation 



\x a b\ Z = X! X! l Xafc ^ ^ 5 ( £a - e b-tv^) S(e a - e F ) 

a b^a 



A dt 



Sx 2 



l(x(0)-x(t)) 2 



cos tot (9) 



where the last line should be averaged over all states 



at the Fermi energy; (5a 



R is the position un- 



(x(0)-x(t)Y 



may be 



certainty in the state a; and 

approximated by averaging over the appropriate classical 
trajectories. 

For a dot in the diffusive regime, where v f t <C R, we 



have t r = R 2 /2D. Then, (x(0) - x(t)) 

as long as t <C r; it grows as 2Dt for r 
finally approaches 2 (5x 2 ) aa for t > tr. Thus 

~ c 2 -%k(u,T R ) 2 for 



L0 < T 



,2_DA 



for 



< tO < T 



It 
-1 



(10) 



where c is a constant which depends on the dot's shape. 
For a roughly circular dot of radius ii,c~ O(l), and we 

ignore it below. The value of |(w x ) Qb | 2 falls off according 
to m) as to increases further. 



2 



For a ballistic chaotic dot, the time scales r and 
tr = R/vf coincide, and D rj vpR/2. The second line of 

Eq. ( |l0| ) does not apply. The maximum value of | Cu^c ) afo | 2 
is RJ VfRA/ttH, obtained when lutr rj 1. From these re- 
sults, we may calculate e so and X 2 (E Z ) = 7 2 |( , y a ;) a h| 2 /A 2 . 
For both types of dots the confinement leads to a S|| de- 
pendence of \ 2 ((jj) for Eztji <C h (see the inset of Fig. 1). 
The maximum value of A 2 is A 2 rj 7 

max nhl\ 

If 7 and 77 are both non-zero, A depends on the direc- 
tion of B\\ within the x — y plane. It is different, e.g., for 
.B||||(110) and £?||||(110), even for a dot which is roughly 
circular. However, the average of A 2 over all directions 

of B\\ will be ( 7 2 + r, 2 )\(v x ) ab \ 2 /A 2 . 

In order to make a comparison to the experiments of 
Ref. M, we consider a ballistic dot where 7,77 are such 
that A(A) <C A c and A max is greater than the crossover 
value A c (see Fig. I). At low temperatures T <§; H/tttr, 
the variance var(<?) should decrease in two stages as the 
Zeeman energy Ez is increased. In the first it would drop 
from 4C% to 2C% over the range < E z < n max(T, HT), 
due to the removal of spin degeneracy of the levels. Then 
var(g) would drop by an additional factor of approxi- 
mately 2, resulting from the turning on of spin-orbit cou- 
pling, over the larger range 7rmax(T, HT) < Ez < Ti/tr. 
(For still larger value of Ez , the conductance fluctuations 
would increase again.) In Ref. [Q, T was comparable to 
h/ftTn. Under such conditions we expect the factor-of- 
four decrease in var(g) to occur smoothly over the range 
< Ez < fi/rR. (Recall that time-reversal invariance is 
broken by B± in all cases). 

A quantitative comparison of this scenario to the ex- 
periment of Ref. Q requires information regarding the 
strength of spin orbit coupling, which we parametrize by 
the dimensionless parameter Q so = (j 2 + t^Y'^Vf / Ep- 
In terms of Q so , for a ballistic dot A max rj Q so Ne^ 4 /3, 
where N e is the number of electrons in the dot. For our 
scenario to be consistent with the experiment, we need 

1 12 

A max to be at least comparable to A c ~ 0.2N c g for the 
large measured dot (where N e — 16,000 and N e g RJ 6) 
and smaller than A c for the small measured dot (where 
N e = 2,000 and N c g rj 6). These requirements suggest 
5 • 10~ 3 > Qso > 10~ 3 . There are additional uncertain- 
ties, however, because our application to ballistic chaotic 
systems of formulas derived for diffusive systems (e.g., 
Eq. |Io|) involved several unknown numbers of order unity. 

Although 7 and rj have been measured previously in 
other GaAs heterostructures, the parameters depend on 
details of the structure, and are difficult to extrapolate 
from one system to another. Values of Q so extracted from 
existing data on GaAs 2DEGs include Q so RJ 1.6 x 10~ 2 , 
from optical measurements in a sample with n = 4x 
10 11 cm~ 2 , and Q so « 5 x 10~ 3 , from Shubnikov-de Hass 
measurements Q, in a sample with n = 1.2 x 10 12 cm~ 2 . 
Magneto-resistance measurements JTu] in 2DEGs extract 



the spin-orbit scattering rate by studying the crossover 
from weak localization to weak anti-localization as the 
density is increased. At the densities where the crossover 
occurs, typically around n ~ 6 x 10 11 cm -2 , values around 
t^ 1 ~4x 10 10 sec _1 are found in the moderate mobility 
samples (I ~ 0.5/um), corresponding to Q so ~4x 10~ 3 . 
Our estimates of Q so for the samples of Ref. pL which 
had n = 2x 10 11 cm -2 , are not incompatible with the 
range of previous measurements. 

The suppression of spin-orbit matrix elements by the 
confinement to a dot affects also the scattering rate due 
to spin-orbit coupling. Note that spin-orbit scattering 
processes do not necessaruly result in a spin flip of the 
electron. The probability of a spin-flip in a spin-orbit 
scattering process depends on the ratio 7/77 and on the 
intial direction of the spin. We focus on the case 77 = 
and initial spin state in the x — y plane, in which half 
of spin-orbit scattering processes involve a spin- flip. We 
also set Ez = for this part of the discussion. The rate of 
spin-flip due to a spin-orbit scattering process of a state 
a) is — = Im S(a, e a ), where E(a, e a ) is the on-shell self 
energy of the state \a) due to spin-orbit scattering events, 
irrespective of the final spin-state. 

To second order in the spin-orbit interaction, the en- 
semble averaged imaginary part of the self energy de- 
pends on I a) through e a . Since the total ImE is twice 
the contribution from the spin flip process, we have 



ImS(o,e) = 2tt^2 |(#so) aT M S (e - e b 

b 



fill 



We neglect the real part of the self energy, which does 
not significantly affect our considerations. 

Due to the finite escape rate T from the dot, which is at 
least comparable to the level spacing A, the ^-functions 
in ( pd| ) are broadened enough to allow the sum to be 
replaced by an integral. Then, in view ofEqs. (|) and 
(|Io|), the on-shell self energy T,(a,e a ) vanishes. We go 
beyond this order, to a self-consistent self-energy, where 



E(a, £ a ) 



deb 



(Hso) a 



A e a - e b - E(6,e a ) - iT 



(12) 



and approximate the solution of (]lj) by substituting ( |Tl| ) 
in its right hand side. In the diffusive limit we find, 



j_ J| (4^ + r) for i| < 1; Tt r « 1 



for ^>1 



(13) 



where — A^D/ti 2 is the spin-flip scattering rate in an 
open system. As expected, when tr 3> the confine- 
ment of electrons to the dot does not significantly affect 
spin-orbit scattering rate. In contrast, for a small dot 
spin-orbit scattering rate is suppressed. When T <C 7^ 
it becomes of fourth order in the coupling constants, 
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^- ~ i7r (JLyi ■ At this order, the smallness of spin-orbit 
matrix elements at close energies is overcome by virtual 
transitions of high energy difference (~ t^ 1 ). 

The suppression of r" 1 in small dots has implications 
for electronic transport through the dots. Normally, for 
an open system, when t^ 1 gets larger than r^ 1 , weak 
localization turns into weak anti-localization, and the 
magneto-resistance becomes positive. The correspond- 
ing criterion for a quantum dot compares the reduced 
r^ 1 with (t^ 1 + T). In fact, the spin-orbit scattering 
rate relevant for transport may be even smaller than the 
already reduced r" 1 given by ([l3]). Even in open 2D elec- 
tron structures in GaAs, there are subtleties due to the 
fact that spin rotations induced by -ff so are correlated 
with spatial displacements of the electron. As a result 
the spin-orbit relaxation rate relevant for weak antilocal- 
ization and for conductance fluctuations at Ez = can 
be smaller than that obtained from the above analyses. 
In fact, when rj — ±7, there is actually no weak antilo- 
calization, if the cubic term, oc v 3 in _ff so , is ignored p2|. 
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FIG. 1. Variance of conductance for the case 

N = 1, t^ 1 — 0, as a function of the coupling A, for T = 
and T = A. Asymptotes show known results for A = and 
A = 00 at T — 0. Data for T = A were scaled by a factor of 
5.5. The inset shows schematic behavior of A as a function of 
the Zeeman energy E z , following Eqs. (|) and @. 

Before concluding, we explain the random matrix cal- 
culations leading to Fig. [1]. These calculations are 
aimed at studying the A-dependence of var(g) in the pres- 
ence of a strong Bu, The Hamiltonian of the closed dot 
was modeled by a 2M x 2M random matrix H of the form 
Hij = rjijKij where 1 < i < M labels states with spin up, 
M + 1 < i < 2M labels states with spin down, rjij = 1 if i 
and j are states with the same spin, rjij = e so — nX/M 1 / 2 
if i and j are states with opposite spin, and K is a random 
Hermitian matrix from the GUE, with the distribution 
P(K) cx e~3 tr K * K . This matches the definitions in the 
text because e so is the RMS value of the matrix clement 
connecting two states in the different spin blocks, and 
A = -k/M 1 / 2 is the average level spacing at the center of 



the band for one block. For the case N — 1, we connect 
the system to "leads" with perfect conducting channels 
at 2 sites for spin up and 2 sites for spin down, and cal- 
culate the 2x2 transmission matrix t for energies near 
the center of the band pT| . 

At T = one finds the conductance by using the Lan- 
dauer formula g — tr tt^ . To obtain results at finite T, for 
each realization of the random matrix we first calculate 
the transmission matrix and thus the T = conductiv- 
ity g (0) {E) for a range of energies E. This conductance 
is weighted by the derivative of the Fermi function and 
integrated to give g(T) = ^g^(E)dE. 

The conductance fluctuations as a function of A for zero 
temperature and temperature A are shown in Figure 1. 
The data for T = A were accumulated from 5000 real- 
izations with M = 60, while the data for T = were ob- 
tained from 10 6 realizations with M = 20. More limited 
calculations at T — with M = 60 showed differences of 
less than 10 % from M = 20. 

It is clear from the figure that the T = A results have 
the same depedence on A as the T = results, and we 
have evidence that this remains true for higher T. The 
reduction in var(g) by a factor w 6 is consistent the theo- 
retical expectation @ that C N {T) w C N {T = 0)7r^r/6T, 
for T > 7rfiF, if dephasing is absent. The middle of the 
crossover occurs roughly at A = 0.1, which is somewhat 
less than the value 0.23 given by Eq. (||). 

Calculations for 2 < N < 4 at T = (not shown) are 
consistent with a crossover value A c scaling as A 1 / 2 , as 
predicted by (||). Dephasing can also be included by us- 
ing a "third-lead" model of the type discussed by Brouwer 

and Beenakker [R. Results for N — 1 and T — seem 

. 1/2 
to show a variation of A c somewhat slower than N L , at 

least in the range 1 < N c g < 4. 

In summary, we presented a theory by which the ef- 
fect of spin-orbit on conductance fluctuations in a quan- 
tum dot depends strongly on an applied parallel magnetic 
field. This theory may well explain the experimental ob- 
servation of Ref. 0] . 
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